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Abstract 

We use a discrete-time formulation to study the asymmetric avalanche process [Phys. Rev. Lett. 
87, 084301 (2001)] on a finite ring and obtain an exact expression for the average avalanche size 
of particles as a function of toppling probabilities depending on parameters \i and a. By mapping 
the model below and above the critical line onto driven interface problems, we show how different 
regimes of avalanches may lead to different types of critical interface behavior characterized by 
either annealed or quenched disorders and obtain exactly the related critical exponents which 
violate a well-known scaling relation when a^2. 
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The physics of driven interfaces has been attracting great attention for many years due 
to its connection with a variety of phenomena such as fluid flow through porous media, 
motion of charge density waves and flux lines in superconductors, etc. 

HQS. The 

results obtained can be roughly splitted into two groups. The first one incorporates growth 
phenomena where thermal or annealed fluctuations affect the dynamics. The most famous 
example is the Kardar-Parizi-Zhang (KPZ) universality class characterized by the roughness 
exponent x — 1/2 and the dynamical exponent z = 3/2 in 1 + 1 dimensions |^|. In the second 
group, an interface moves under the action of an external driving force F, which competes 
with pinning forces due to quenched medium inhomogeneities. At some critical force F c the 
transition occurs from the totally pinned state to the state where the interface moves with 
a velocity v, which shows a power law decay v ~ (F — F c ) when F approaches the critical 
point F c from above [5( . 

The scaling properties of interfaces in the quenched case strongly depend on the medium 
isotropy. The isotropic medium produces a rough interface, which is believed to obey the 
quenched KPZ equation with a nonlinear term vanishing at the depinning threshold, F — > F c 
ja, Q]. In the anisotropic medium, the quenched KPZ equation holds only for a definite ori- 
entation (hard direction) of the interface, yielding a divergent nonlinear term. The tilt from 
the hard direction generates the gradient term in the equation which breaks the translation 
invariance. The existence of such a term gives rise to another universality class, tilted in- 
terface (TP class , characterized in 1+1 dimensions by the exact exponents \ = 1/2, z = 1, 
and 9 = 1 0, li- 
lt has been noted that thermal fluctuations below the depinning threshold can initiate 
avalanches, which lead to infinitely slow creep of the interface jjj- Far below depinning 
threshold, the avalanches lead to local advances of finite interface segments, so that they 
can be considered under coarsening as a local annealed disorder in an appropriate time scale. 
Approaching the depinning threshold, the local avalanches increase up to the system size 
and contribute to a global interface depinning which is controlled by the quenched disorder. 
It is the aim of this Letter to treat the effects of annealed and quenched disorder in a frame 
of a unified model. 

Due to the well known correspondence between growth processes in (1 + 1) dimensions 
and one-dimensional lattice gases ^J, we can interpret the avalanche dynamics of interfaces 
in terms of avalanches of particles similar to those appearing in one-dimensional random 
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sandpile models jo]. To this aim, we can use the Asymmetric Avalanche Process (ASAP) 
which has been recently proposed and solved by the Bethe ansatz method ^(|. Having 
formal origin in the asymmetric exclusion process h|, the ASAP has completely different 
dynamical behavior specific for systems with avalanche dynamics which has been widely 
investigated in the context of self-organized criticality 1^ . 

The ASAP is formulated as a system of particles diffusing on the lattice. The interac- 
tion makes particles to accelerate if they are in the same site of the lattice, giving rise to 
avalanche-like relaxations of domains with higher density of particles. Such a dynamics de- 
pends on transition rates and may lead to different regimes of the particle flow, intermittent 
or continuous. In the former case, fast avalanche processes alternate with long periods of 
relatively slow diffusive motion and the average avalanche size of particles is finite being 
measured in the diffusive time scale. This average avalanche size di verg es in the thermo- 



dynamics limit if the density of particles approaches a critical point [10( • The divergency 
means that the diffusive time-scale is no more appropriate for the description of the system, 
because the avalanche processes become dominating. Internal avalanche time should be used 
in this case to quantify the dynamics in the system. 

The study of the intermittent flow phase by the Bethe ansatz on a lattice of finite size 
N Q| shows that the total distance traveled by particles has the long time asymptotic 
(Y 2 ) c ~ tN 3 / 2 , which is equivalent to finding the KPZ dynamical exponent z = 3/2 (jj, \$. 
. t will be clear below that in the phase of continuous flow, Maslov-Zhang (MZ) arguments 
^| based on the mapping of the system onto the random walk problem show that the 
model belongs to the TI universality class at least in the vicinity of the transition point. In 
this Letter, we use the equivalence between the driven interface problem and the ASAP in 
discrete time to study the crossover between these two types of disorder. Considering the 
ASAP on the lattice, we obtain the exact expression for the average avalanche size in two 
different regimes and describe scaling properties of the system at the transition point. 

Consider the system of p particles on a one-dimensional lattice of N sites with periodical 
boundary conditions, i.e. a ring. The system evolves by discrete time steps At < 1 according 
to the following rules. I. If all particles occupy different sites of the lattice, then one 
randomly chosen particle moves one step to the right neighboring site with probability At 
for one time step. II. If there are n > 1 particles at the same site, referred to as an active 
site, the diffusive motion is replaced by the avalanche dynamics: either all n particles move 
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to the next right site at the next time step with probability \x n or n — 1 particles move one 
step to the next right site with probability 1 — fi n , while one particle remains unmoved. 
Thus, we use a random sequential update and separate in time the avalanche and single 
particle dynamics. The continuous time limit, iV — > oo and At — > 0, coincides with the case 
considered in ^(|, as only the terms of order At survive in the master equation. 

To map the system to the interface problem, we associate the occupation number n» 
on a site i with an interface height increase rij = h (i + 1) — h (i) so that the number of 
particles determines the tilt of the interface. To combine periodic boundary conditions 
with the interface tilt we imply helicoidal boundary condition to the interface by putting 
h (i + N) = h (i) + p. 

A state of the system at time t is given by probabilities of all possible configurations of 
particles on the lattice. Such probabilities obey the master equation 

p t+ At(c) = j2nc,c)p t (C) (i) 

{C} 

where T(C, C) denotes the transition probability from the configuration C to C and the 
summation is over all possible configurations of particles. Despite the ambiguous explicit 
form of Eq.([T|). the form of its stationary solution is rather simple and can be obtained by 
direct solution under an assumption of stationarity. In the stationary state, the probability 
of a configuration depends only on the number of particles in the active site rather than on 
occupation of the other sites of the lattice. Let be the probability of any configuration 
with n particles at the active site. Then, the stationary solution of the system (JIJ) is given 
by the following recursion relations 

7> [n+1] =7> [n W(l-/Wi), (2) 

where n — 1, 2, . . . , p — 1, and fii = At /p. The only undetermined constant here is the 
probability of stable configurations with no active site, which is defined from the normaliza- 
tion condition X]{c} = 1. 

As it follows from the Eq.(j2J, when n particles leave an active site at some step of an 
avalanche, all configurations of the remaining p — n particles at the lattice occur equally 
likely. This form of the stationary measure is an analogue of the equiprobability of all 



configurations in the exclusion process 



in Ref. and consider the expected number (n) of events when n particles leave an 



111 ]. Then, we can proceed in a way similar to that 



active site at the k-th step of an avalanche. The number of particles, n, outgoing from an 
active site at every time step performs Markov random walk 



Pk+i(n) = ^ftH W m,n (3) 



with transition probabilities 



n — 1 
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for n > 1 and w M = (p - ^) (1 - p 2 ), ^2,i = (l - (p - f )) (1 - P2), where p = p/N. 
Equation (j^J generalizes corresponding equation in [10] explicitly taking into account the 
effect of the finite size on the transition probabilities. This enables us to consider the ASAP 
both in the intermittent and continuous flow regimes. Instead of Pk{n), we can look for 
the total expected number of these events during the whole avalanche P{n) = J2T=i Pk{n), 
which satisfies much simpler equation P(n) = S Ui i + ^2 m w mn P(m), supplemented by the 
initial condition P (1) (l — ^jf) = 1. The solution of this equation can be found from the 
detailed balance condition, P in) u> n ,n+i — P (n + 1) w n+ i jTl . The result is 



P(n) 



NT(p)T(N -p+l) A Pi 



r(p - n + 1)T(N - p + n + 1) -fJL 1 _ ^ 



for n > 1. The expected total number V of particles spilled during an avalanche, i.e. the 
average avalanche size, is given by 

p 

V = J2 nP ( n )- ( 5 ) 

n=l 

As a single-particle jump which can trigger an avalanche occurs with probability At, the 
average number of steps of particles for time At is AtV, if one neglects the time taken 
by the avalanche itself in comparison with the average time interval between avalanches. 
Therefore, if At — > and V is finite, we may consider V as an average velocity of particles 
or the velocity of the associated interface in the continuous diffusive time scale. 

Consider the asymptotic behavior of the average avalanche size, V, in the limit iV — > 
00, p — » 00, p/N = p = const. To this end, we have to specify the toppling probabilities fij. 



We suppose that \ij tends to a constant limit /i^ when j — > oo, and define the function / (n) 

through the relation /i„/(l — u n ) = e^ n ^Poo/(^ — Poo), assuming that it has the expansion 

/ (n) ~ (a — 2)/n + 0(l/n 2 ). It turns out that a is the only parameter which is responsible 

for the subcritical singularity of the average avalanche size. Below we imply a to vary in 

the range 1 < a < 3. In the thermodynamic limit we have the leading term of the n/N 

( \ n 1 

expansion V ~ Y^=i n0l ~ l ( (l-^f °°)(i-p) ) > which converges for a > if the density of 
particle is less than the critical value p < p c , where p c = 1 — p^. Close to the critical point, 
V diverges with the critical exponent a, 

V~{ Pc - P y a . (6) 

Above the critical line, we use the Stirling formula to approximate gamma functions in 
Eq. (@J). The leading order of the expansion shows that the function P (n) has a maximum 
at the point 

n = N (p - p c ) . (7) 

The formula (J7|) has a transparent physical meaning. Indeed, n is the most typical value 
of the number of particles at an active site during an avalanche, which ensures the density 
of particles at the other sites of the lattice to be equal to p c . In other words, the avalanche 
dynamics in the system above the critical point is self-organizing in a sense that it maintains 
parameters of the medium outside the active site at the critical point. In the vicinity of no, 
the function P(n) shows a Gaussian-like form with the power law prefactor n a ~ 2 . Replacing 
the sum (jSJ) by the integral, which can be evaluated in the saddle point approximation, we 
obtain the leading part of V, exponentially growing with iV for p > p c 

V ~ (p - p c ) a ^ N a ^/ 2 a N (8) 

where a = (p/p c ) p ((1 — p)/{\ — p c )) 1 ~ p . Thus, the average avalanche size either exponen- 
tially diverges with iV for p above p c or tends to a finite limit below p c . These two regimes 
are separated by the point p = p c , where the power law growth of the avalanche size takes 
place 

V c ~ N a/2 . (9) 

The crossover from sub- to super-critical regime through the critical point can be viewed 
in the function obtained from another expansion of Eq. (j3J) for 1 n N, namely, 
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V ~ V c g a {t), where 



2e t2 



r (a/2) 



9* (t) = T^-nr, I dx(x + t) a - 1 e~ x2 (10) 



is a function depending on the parameter a and variable t = ^ N (1 — p) p/2 In ^Zpff i which 
shows a distance from the critical point. Near critical point t ~ (p — p c ) y/~N. The function 
g a (t) is equal to 1 when t = 0, and decays as \t\~ a when t — > — oo, thus eliminating the 
dependence of V on N in the subcritical regime. Next, we consider particular examples of 



toppling probabilities discussed in 



In the integrable version of the ASAP, the toppling probabilities are defined by the re- 
cursion relations p 2 = Pn+i = Ml ~ /V)- m ^ ne notations used, this case corresponds to 
/^oo = /J-/ (1 + /x) and a = 2. The sum in © can be reexpressed in terms of the infinite sum 
of hypergeometric functions 2F1 (a, b, c, d) 



v =^tf5rrE^ -(-M fc+1 x (11) 



li N-p+1 



\k+l 



2 F 1 (2,l-p,2 + iV-p,(-/i) 



k+l 



This form is useful because all terms of the sum are of the same structure. The terms with 
even k look like the average avalanche size of the ASAP with particular toppling probabilities, 
p n = p k+1 , independent on n, hence having its own critical density p^ — 1/ (l — (— /i) fe+1 j, 
while the odd k-th terms are not singular in all phase space. Therefore, for pi 2ra ' < p < p^ n+2 ^ 
the growing part of the average avalanche size is given by the sum of n + 1 exponentially 
growing summands similar to one in Eq. (JSJ) up to the change of p c to p£ k \ (k = 0, 2, ... , 2n) 
and a = 2. In the thermodynamic limit t 
point p c = 1/ (1 + p) , the exact equality 



re term with k = is dominating. At the critical 
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2 F 1 (2, 1 - p, 2 + pp, -p) = (1 + + /i), (12) 

leads to the surprising result,^ = iV + 0(1), which shows that the growing part of the 
average avalanche size does not depend on p in the integrable case. The expansion below 
the critical point, p < p c , confirms the thermodynamics value and 1/N correction obtained 
from the Bethe ansatz solution 

HQ. 

For two-parametric ASAP [l0( with toppling probabilities defined by p n /{l — p n ) = 
2p{n — 2 + a)/((l — p)na), we have p^ = 2pj (a + 2p — ap) and 

^ = F^TT 2F '( a • 1 ~ ^, • 2 + iV ~ ^, '^))• (13) 



When a = 2, we have p n = p^ = p, p c = 1 — /i, and V c = p. 

Let us consider a connection of these results with MZ arguments jjj. Considering the 
form of transition probabilities of Eq.fjHJ) for p = p c , one can see that the height of an active 
site in the ASAP also performs the simple random walk when a = 2. However, for a / 2 
the normalized transition probabilities, defined as p(n — > nil) = w n , n ±i/ \w n ,n+i J r w n,n-i)i 
contain the nonuniform bias term, p(n — > n i 1) = 1/2 l±(a-2)/ (2n) + 0(l/n 2 )), which 
corresponds to the process first considered by Gillis |l6|. Depending on whether a is larger 
or less than two, it can be either positive or negative, respectively, enhancing or suppressing 
the avalanche spreading. 

As usual, we assume the following scaling ansatz for avalanche size and time distribution 
P(s) = s~ T g(s/so) and P(t) = t~ Tt g(t/t Q ) in the vicinity of the critical point, where g(x) 
is a scaling function, and so,to are time and size cutoffs. The time cutoff to also plays 
the role of correlation length. Knowing asymptotics for mean conditional time of the first 
return for the Gillis random walk [17], we can obtain the avalanche time critical exponent, 
Tt = 5/2 — a/2 for 1 < a < 3. To probe the avalanche size statistics we should note that 
unlike the exponent r t , the dimension of avalanches D does not depend on a coinciding 
with unbiased case, which can be directly checked by multiplying the original equation © 
by n 2 and summing by parts. This yields (n 2 ) ~ t and hence D = 3/2 and r = 2 — a/3. 
The other critical exponents characterizing underlying interface dynamics remain unchanged 
comparing to unbiased MZ case. Particularly, the exponent characterizing the correlation 
length below the critical point, defined in our case by to ~ (Pc — p) u , can be uniquely fixed 
as v = 2 by requiring the obtained exponents to be consistent with the average avalanche 
size below p c , Eq.flOJ). In the critical point one gets from the Eq.® to ~ N irrespectively 
of a, i.e. the time cutoff grows as a system size, which determines the dynamical exponent 
2 = 1. We should note that although the correlation length demonstrates universal behavior 
independent on a, both in the subcritical (KPZ) and critical (TI) regime, the results for the 
mean avalanche size in the intermediate regime (fTT)|) show the remaining dependence on a. 
The exponent 6 = 1 just above the depinning threshold follows from the Eq. fTj) . which gives 
a characteristic increase of the interface height after one avalanche passage proportional to 
(p — p c ). The roughness exponent x = 1/2 is due to equiprobability of particle configurations 
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established above. We should note that the case a^2 violates the scaling relation 

r = l + (d-l/u)/D (14) 

valid for a wide class of the models [9] probably because the interaction is not uniform with 
respect to the height of an avalanche. Rigorous explanation of this violation is a matter of 
further analysis. 

In summary, we have found that the asymmetric avalanche process (ASAP) [n| incorpo- 
rates annealed KPZ and quenched TI interface dynamics for p < p c and p > p c , respectively. 
All critical exponents characterizing the interface behavior in both classes are obtained ex- 
actly and shown to coincide with those known before. Nevertheless, the exact calculation 
of the average avalanche size and random walk arguments shows that critical exponents 
of the avalanche distributions are the continuous functions of the parameter a responsible 
for the asymptotics of toppling probabilities and the scaling relation between avalanche- size 
exponent r and correlation-length exponent v 

[9] valid for many other models is violated in 

the ASAP with a ^ 2. 
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